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SUMMARY 


NEWSUMT is a computer program written in FORTRAN subroutine 
form for the solution of linear and nonlinear constrained or un- 
constrained function minimization problems. The basic algorithm 
is the sequence of unconstrained minimizations (Ref. 1) using 
the modified Newton's method (Ref. 2) for unconstrained function 
minimizations . 

Problems must be formulated in the following form: 

Minimize F(X 1/ X_,...X ) 

x z n 

Subject to g q (X r X 2 ,...X n ) 2 0 q = 1,2, ...Q 

The user must provide a main program which calls subroutine 
NEWSUM and also subroutine ANALYS which computes the function 

F (X) and g q (X) . If analytic gradients of these functions 
are available, the gradients should also be computed by ANALYS; 
otherwise the gradients will be computed by finite differences. 
Even if constraint functions or the objective function are not 
defined for certain values of the design variables, artificial 
definitions must be specified so that all functions are defined 
and differentiable over the entire design space. 

This report describes the use of NEWSUMT and defines all 
necessary parameters. Sufficient information is provided so 
that the program can be used without special knowledge of non- 
linear mathematical programming methods. 
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I . INTRODUCTION 

NEWSUMT is a computer program written in a FORTRAN sub- 
routine form for the solution of linear and nonlinear inequality 
constrained or unconstrained function minimization problems. 

The purpose of NEWSUMT is to determine the values of a set of 
variables X (a vector of real variables, X x ,X 2 , . . . X^^.) that 
minimize a function F(X) subject to a set of inequality con- 
straints g (X) i 0, q = 1,2,...,Q. NEWSUMT was originally 
developed as an optimizer for sizing minimum weight finite 
element structural systems in the ACCESS-2 computer program 
(Ref. 3). However, it is a general purpose optimizer that can 
be used for solving a wide variety of numerical optimization 
problems. It treats inequality constraints in a way that is 
especially well suited to engineering design applications. In 
a structural design application, for example, X could represent 
dimensions of structural members, F(X) could be the weight of 

the structure, and the g (X), q = 1,2,...Q, could express stress. 

Si 

buckling, or other types of behavioral constraints. 

When using NEWSUMT the optimization problem must be formu- 
lated in the following form: 

Minimize the objective function F()t) 

Subject to inequality constraints 

g q (X) k 0, q = 1,2, ...Q ( 1 ) 

Xj (L) < Xj £ Xj (U) , j = 1,2,... NDV 
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where the functions F(X) and g (X) are continuous and differ- 

Si 

entiable real functions with respect to the design variables 

x j > j = 1/2,...,NDV. The user must supply a main program which 

calls NEWSUMT and also a subroutine which is called by the 

NEWSUMT program to evaluate functions F(X), g (X) and, if avail- 

able, the derivatives of F and/or g with respect to the vari- 

ables , X j . The user specifies an initial design by assigning 

certain numerical values to (X 1 ,X 2 , . . . ,X NDV ) ; the NEWSUMT program 

then systematically modifies these values generating a sequence 

of vectors X such that F(X) decreases and none of the inequality 

constraints are critical. This sequence of vectors X converges 

to a solution X where all the inequality constraints are satis- 
• ->■ * 

fied and F(X ) is at least a local minimum. 
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II. 


MINIMIZATION ALGORITHM 


In this section, the algorithms used in the NEWSUMT code 
are explained. All of these computations are performed inter- 
nally. They are described here for completeness. 

The minimization algorithm used in NEWSUMT is a sequence 
of unconstrained minimizations technique (SUMT) , Ref. 1. Major 
features of NEWSUMT which distinguish it from the original 
formulation (see Ref. 1) include: 

(a) A modified Newton's method is used in the direction 
finding part of the unconstrained minimization- In 
this method, second derivatives of the constraints 
are approximated by expressions involving only the 
first derivatives (Ref. 2) 

(b) An extended interior penalty function formulation. 
This type of penalty function combines the features 
of interior and exterior penalty functions. That 
is, although initial designs may violate the con- 
straints, subsequent designs satisfy the constraints 
and tend to be noncritical (Refs. 4,5). The trans- 
sition point control parameter, which is a critical 
factor for numerical stability when using extended 
penalty functions, is selected so that the one 
dimensional search problems generated usually have 
their minimum points inside the feasible region 
(Ref. 6). 

Since the detailed mathematical aspects of the SUMT 
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algorithm and its variations are beyond the scope of this report, 
only matters of critical importance to understanding the funda- 
mental procedure used in the NEWSUMT program will be described 
in the following subsections. 

2 . 1 SUMT Approach 


The SUMT algorithm transforms the inequality constrained 
problem defined by equations (1) into a sequence of unconstrained 
problems. To accomplish this transformation, a compound func- 
tion <p(X,r ) is introduced. The compound function used in 
NEWSUMT is defined as 



( 2 ) 

In the transformed problem, <f>(X,r ) is minimized with respect 

P 

-)*• _ 

to X for a sequence of decreasing values of r , which is called 

P 

the penalty multiplier. Because r is being decreased, the 

P 

contribution of the penalty function is being reduced and the 
solution to the transformed problem is converging toward the 
solution of the initial problem defined by Equation (1) . Fun- 
damental theory, convergence characteristics and various methods 
for solving SUMT-type problems are discussed in Ref. 1. 

Optimization technology developed during recent years 
has greatly improved the computational efficiency and numerical 
stability of SUMT-type formulations. NEWSUMT incorporates 
many of these new features and they enhance its performance. 
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2 . 2 Extended Penalty Function 


It must be noted that the composite function given by 

Eq. (2) is only defined for that portion of the design space 

where all the inequality constraints are satisfied. This is 

extremely inconvenient because it requires that an initial 

design which satisfies all constraints with a reasonable margin 

be available. Furthermore, since cf>(X,r ) is undefined in the 

P 

infeasible region, the one dimensional minimization algorithm 
becomes complicated and inefficient. Using the extended pen- 
alty function concept, Eq. (2) is modified as follows 

+ + T Q NDV 1 

<J)(X,r ) = F(X) + r l H (X) + l [L (X . ) + U(X )] (3) 

P P Lq=l q j=l ^ 3 J 

where 



3g q (x) 

e 


+ 3 


g q (x)^ e 


(4) 


g q (XX 


E 


U(X j ) 



3 (X.-X . (L) ) 
2 L_ 


+ 3 


X j ^=- (unbounded) 


X . 


3 


(L) 

j 


(5) 


X. 


3 


< 


X. (L > 

3 
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and e denotes the transition parameter ( Sec. 2.3). The advantage 

of the extended penalty function concept becomes apparent upon 

examining Fig. 3. The original penalty function defined by 

Eq. (2) is shown as a broken line which approaches infinity as 

g (X) approaches zero from the positive side and it is not de- 

fined in the infeasible region where g (X) <0. On the other 

hand, the penalty function H (X) given by Eq. (4) is defined for 

negative values of g (X) and it is a smooth function (continuous 

si 

up to the second derivatives) of g (X) . As a consequence, 

g 

Hg(X) is a well behaved smooth function within the region where 

g (X) is defined and it is also a smooth function of (X). An 
Si. 

interesting interpretation of quadratic extended penalty function 
is given in Appendix A, where it is viewed as a combined interior- 
exterior penalty function approach. 

2 . 3 Transition Parameter 

The transition parameter (e) , introduced when defining ex- 
tended penalty functions (Eqs, 4 through 6) , is selected 
initially by the user (default value is 0.1) and may be given 
as input data. As soon as a design which satisfies all con- 
straints is found after two or more unconstrained minimizations, 
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the value assigned to e is automatically estimated by the 
method set forth in Ref. 6. This method guarantees that the 
transition parameter is chosen as large as possible (to main- 
tain numerical stability) while at the same time ensuring the 
minimum remains inside the feasible region. Once e is deter- 
mined automatically, the coefficient C which relates e with r 

Jr 


E = c/F" (7) 

p 

is computed; then Eq. (7) is used thereafter to compute e from 
the response factor r^. 

2.4 Modified Newton's Method 


Minimization of $ (X, r ) with respect to X involves repeated 

P 

application of two basic steps; 

Step 1. Find a direction vector S along which the design is 

-V 

modified starting from the current design Xq . A new 
design X is given by 


X = X Q + aS (8) 

where a is a scalar variable that governs the move 
distance in the design space. 

Step 2. Find the value for a so that the composite function 

<f) (X, r ) is minimized along the direction S. 

P 

The modified Newton's method discussed in this subsection deals 

“ 4 " 

with the procedure used to find a direction vector S (Step 1) . 

-4 , , , 

In Newton's method, the direction vector S is given by 
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S = - [J]" 1 V<}>/ ||[ J] ” 1 V<f> || (9] 

where [J] is a NDV x NDV matrix with (i,j) element defined by 

J ij = 3X. 8X . (X ' r p* ^ 1( 

J ID F 

Equation (10) may be evaluated by differentiating Eq. (3) . 


2 

3 H 

2 _ ( h - 

3X.3X. (A> 

i D 


2 P g q 9g q _ q 

3X. ax. "q 3X 

i i ^ 


!f2_l 

. ax . 

i i J 


g q (x) > e 


rag ag a 2 g q . 

! Lax q ax' 3 + (2g q " 3£) ax. ax ~V Z 
L i d ^ i i J 


g q (x) < e 


Following the suggestion in Ref. 2, h (X) is simplified by 

2 

neglecting the terms involving 3 g^ (X) /3X^3X ^ , hence 


ax. ax . 

i D 


(X) = 


3g 3g _ 

2 _Ia /q 3 

ax. ax. /y q 

ID 

3g 8g 0 

2 —3. _q. /e 3 

ax. ax. 7 

1 D 


g q (x) 2 e 


g q (X) < e 


This approximation is justified qualitatively in the following 
manner. For critical constraints, g (X) is small, therefore 

q 


3g 3g 8 g 

2 — 3 . _ 13 . >> n 13 
ax. ax. g q ax. ax. 

1 D M i D 


assuming that g (X) is a smooth function of X. For noncritical 

constraints, g (X) is large, thus the entire term corresponding 
Si 

3 

to such constraint is small (due to g^ in the denominator) 
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compared to those associated with critical or nearly critical 
constraints . 

The modified Newton's method uses these approximate con- 
tributions to the Hessian matrix (see Eq. 10) in computing the 
direction vector S. For optimization problems involving a 
large number of complicated nonlinear constraints, experience 
confirms that this approach is efficient and generates good 
quality direction vectors so that only 4-6 one dimensional 
minimizations are sufficient for each unconstrained minimization, 
regardless of the number of design variables. If it is observed 
that the direction vector S found by the modified Newton's 
method does not decrease the composite function <f> , S is replaced 
by the direction of the steepest descend, i.e. - V<J>/ ||V<|) ||. This 
is experienced occasionally due to numerical ill-conditioning 
of the approximated Hessian matrix. 

2 . 5 One Dimensional Minimization 

As mentioned in the previous subsection, it is necessary 
to find the value of a for Eq. (8) such that <f>(X,r ) is mini- 
mized along the direction S. This is achieved by first trapping 
a minimum in a finite interval and subsequently by applying 
the golden section algorithm to determine a . with sufficient 
precision. First the move distance implied by the modified 
Newton method is determined. Then the distances to the hyper- 
planes representing linear constraints are calculated. The 
first trial stepsize is then taken as the smallest of the fore- 
going move distances. If the first trial design gives a smaller 
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value for (p than the initial design, then the step size is 
increased by 2.6180 but it is not allowed to exceed the dis- 
tance to the nearest linear constraint hyperplane. This process 

i 

is continued until <J> becomes greater than the previous evalua- 
tion, thus a minimum is trapped within a finite ct interval. 

For the example shown in Fig. 5, a minimum is trapped after 
the third trial design is evaluated. Then the golden section 
algorithm is activated to find a minimum between a„ and a_. 

^ J 

2 . 6 Convergence Criteria 

The golden section algorithm used to calculate a in one- 
dimensional minimizations is terminated if the maximum relative 
difference among the four function values used in the current 
step of golden section is smaller than the specified value 
(EPSGSN) . 

An unconstrained minimization is judged to be converged 
if two successive one-dimensional minimizations do not improve 
the compound function more than the specified fraction (EPSODM) . 
Then, r^ is reduced by r^^r^x RPCUT, and another unconstrained 
minimization is initialized. The entire process is judged to 
be converged if two successive unconstrained minimizations do 
not improve the objective function more than the specified 
fraction (EPSRSF) . 
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III. MAKING NEWSUMT OPERATIONAL 


In order to become familiar with the NEWSUMT optimization 
code, the following steps are recommended. 

1) Obtain the source program deck and sample problem decks. 

2) Read Sec. IV. of this manual for example problems. 

3) Execute the two sample problems. Both problems are 
self-contained and need no input data. 

4) Read this entire manual. 

5) Devise several two to five variables unconstrained 

and constrained minimization problems and solve them 
using NEWSUMT. If the correct optima can be deter- 
mined analytically compare these with optima obtained 
using NEWSUMT. • ' 

6) Experiment with the various options available. 
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IV. PROGRAMMER'S GUIDE 


4 . 1 Main Program 

The overall organization is shown schematically in Fig. 1. 
Usually, all data must be read in the main program (or at least 
before NEWSUM is called) . The primary subroutine of the NEWSUMT 
optimizer is activated by: 

CALL NEWSUM* (name 


where 


* 

,BL 

,BU 

, DDOBJ 

,DG 

,DH 

,DOBJ 

* 

, FDCV 

,FMIN 

,G 

, GB 

/Gl 

,G2 ,G3 

* 

, OBJ 

, OB JMIN 

,s 

, SN 

,x 

,X0 

* 

, IIK 

, ILIN 

, ISIDE 

,NI 

,N2 

,N3 ,N4 

* 

, RAN 

, NRANDM 

, IAN 

,NIANDM) 



name 

z 

Name of 

main subroutine 

of the analys 




program 

which is 

called 

by NEWSUMT 


BL(N1 ) 
BU(N1 ) 
DDOBJ (N3 ) 


DG(N4) 


program. 6 characteris or less. 

: Lower bounds imposed on design variables. 

: Upper bounds imposed on design variables. 

: Second order derivative matrix of the 
objective function. Only upper triangle 
is stored in a vector form; from top to 
the diagonal for each column. 

: Gradient of constraints . Note that the 
storage scheme is : 


g l , l g 2 , 1 ' * * * g NTCE , l g l , 2 * * * g NTCE , 


NDV 


*Variables with single underline should be assigned or initialized 
in the main program. Variables with double underlines are eval- 
uated in the analysis program. 
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where 


DH(N1) 


DOBJ(Nl ) 

FDCV(Nl) 


FMIN 
G (N2 ) 

GB (n2) 

G1 (N2) 

G2 (N2) 

G3 (N2) 

S (Nl) 

X(N1) 

X0 (Nl ) 
IIK (Nl ) 
ILIN (N2) 


ISIDE(Nl) 


q/i 


3g 

= _£ 
ax. 
1 


Internally necessary array (gradient of 
the compound function) . 

Gradient of the objective function 
Stepsize of the forward finite difference 
steps for each design variable. 

Minimum of the composite function. 
Constraint values. 


Internally necessary array to store 
constraint function values. 


: Internally necessary array to store the 

direction vector. 

: Internally necessary array to store 

: alternative design, 

: Initial design variables 

: IIK (i) =i (i+1) /2 

: Linear constraint indicator 

=1 for linear constraints 
=0 for nonlinear constraints 
: Side constraint identification code 

=+l for lower bound only 
= 0 no side constraint 
- =+2 for upper bound only 

=+3 for lower and upper bounds. 
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N1 

• 

• 

=NDV Number of design variables 


N2 

• 

9 

— NCON Number of constraints 




(=1, if NCON=0) 


N3 


NDV (NDV+1) /2 


N4 

• 

9 

NDV*NCON (=1, if NCON=0) 


RAN (NRANDM) 

9 

9 

Real array which may be used in 

the 



analysis program written by the 

user. 

N RAN DM 

9 

9 

Dimension of RAN 


IAN (NIANDM) 

9 

9 

Integer array which may be used 

in the 



analysis program written by the 

user. 

NIANDM 

9 

Dimension of IAN. 



Note that all the variables or arrays with single underline 
must have their values assigned in the main program before the 
NEWSUM subroutine is called. Arrays with the double underlines 
receive specific values in the analysis program written by the 
user. 

In addition to the variables transferred through the 
argument list of the NEWSUM subroutine, a labeled common block 
is used to transfer optimizer control parameters. The common 
block /CONTRL/ must be declared in the main program as 
COMMON/CONTRL/C , EPSGSN,EPSODM,EPSRSF,G0 ,P 
* , RA , RACUT , RAMIN ,STEPMX, 

,IFD , JRPINT , JSIGNG ,LOBJ ,MAXGSN ,MAXODM,MAXRSF 
,MFLAG ,NDV , NTCE 

where C : Not required to assign values. 
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EPSGSN 


EPSODM 


EPSRSF 


G0 


P 


: Convergence criteria of the golden 

section algorithm. For the golden 
section algorithm, there are always 
four function values to be compared 
to each other. If the sum of relative 
difference of these four values is 
smaller than EPSGSN, then the golden 
section is terminated. Default=0. 001. 

: Convergence criteria of the unconstrained 

minimization. If 3 successive one- 
dimensional minimizations do not achieve . 
relative improvements of the composite 
function or more than EPSODM, then uncon- 
strained minimization is judged to have 
converged. Default=0. 001. 

: Convergence criteria for the overall 

process. If three successive uncon- 
strained minimizations do not achieve 
relative improvements of the objective 
function of more than EPSRSF, then the 
NEWSUMT terminates the search process and 
returns control to the main program. 
Default = 0.005. 

: Initial value of the transition parameter. 

Default=0. 01 

: =0.5. Default=0.5. Do not change 
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RA 


RACUT 


RAMIN 


STEPMX 


IFD 


: Penalty multiplier r . Required if MFLAG=1 . 

IT 

Default=l. 0. 

: Penalty multiplier decrease ratio. 

r p+1 =r p x RACUT. Default=0. 1 

: Lower bound of the penalty multiplier. If 

this is zero, numerical instability or 
excessive number of iterations may take 
place. Default=10 13 . 

: Maximum bound imposed on the initial step 

size of the one-dimensional minimization. 
Default=2 . 0. Note that the direction 
vector S is normalized so that § T *^=1.0 
prior to each one-dimensional minimization. 

: Flag for finite difference gradient control. 

= 0 All gradient information must be com- 
puted by the user's analysis program. 

The analysis program should accept 
INFO=3,4, and 5. Default. 

> 0 Use default finite difference step- 
size (0.01) 

< 0 Use user supplied finite difference 
stepsize. 

FDCV(i), i=l, NDV must be specified in the 
main program. 

=1 Gradient of objective function must be 
computed by finite difference. 
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JPRINT 


JSIGNG 

LOBJ 


MAXGSN 


MAXODM 


—2 Gradient of all constraints (including 
linear constraints*) must be obtained 
by finite difference. 

=3 Gradient of nonlinear constraints must 
be by finite difference. 

=4 1 and 2 combined. 

=5 1 and 3 combined 

: Printout control parameter. 

=0 Print initial and final designs only 
=1 Print brief results of analysis for 

initial and final designs together with 
minimal intermediate information. 

(the default option is 1) . 

=2 Detailed printing. 

=3 Debugging printing. 

: Not used. 

: Flag for linear objective function. 

=0 F (X) is nonlinear function of 

(X 1 X 2- * * X NDV ) 

=1 F(X) is a linear function of 
( X 1 X 2 * * * X NDV^ 

: Maximum allowable number of golden section 

iterations. Default=20. 

: Maximum allowable number of one-dimensional 

minimizations per unconstrained minimization. 
Default=6 . 
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MAXRSF : Maximum allowable number of unconstrained 

/ 

minimizations. Default=15. 

MFLAG : Flag for penalty multiplier initialization. 

=0 ; initial r is computed by NEWSUMT. 

IT 

=1 ; r specified by the main program is 

ir 

used as the initial value. 

NDV : Number of design variables. 

NTCE : Number of constraints considered. 

4 . 2 Analysis Program 

The primary subroutine of the analysis program supplied by 
the user should have the following arguments: 

SUBROUTINE name (INFO, X, OBJ, DOBJ, G, GB, DG, Nl, N2 , 

N3 , N4, RAN, NRANDM , IAN, NIANDM) . . 

name : Subroutine name which is identical to the 

first argument of the CALL NEWSUM state- 
ment issued in the main program. 

INFO : Control parameter 

=1 ; evaluates objective function only. 

=2 ; evaluates all constraint functions. 

=3 ; evaluates gradient of objective 
function, 

=4 ; evaluates gradients of nonlinear 
constraint functions. 

=5 ; evaluates gradient of linear con- 
straints only. 

X (Nl ) current design variables. 
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OBJ 

DOBJ(Nl) 
DOBJ (N3 ) 

G(N2) 

GB (N2 ) 
DG(N4) 


objective function value, 
gradient of objective function, 
second derivatives of objective 
function. 

constraint function values, 
usually not used. 

gradient of constraint functions stored 
in a vector form. 

DG(l)=9g 1 /3X 1 

DG(2)=9g 2 /9X 1 


DG(N2)=3g N2 /3X 1 

DG(N2+l)=9g 1 /3X 2 


N1 

N2 

N3 

N4 


DG(N4)=3g N2 /3X N1 
maximum number of design 


variables 
maximum number of 
constraints 
N1 (Nl+l)/2 
N1*N2 


Parameters 
used for 
variable 
I dimensioning 
I of primary 
arrays. 


RAN (NRANDM) 
NRANDM 
IAN (NIANDM) 
NIANDM 


real array allocated by the user, 
dimension of RAN. 

integer array allocated by the user, 
dimension of IAN. 
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Fundamental structure of the analysis program which the 
user must supply is shown in Fig. 2. 


4 . 3 Description of NEWSUMT Subroutines 


NEWSUM 


BLOCK DATA : 


DIRCTN 


FUNCTN 


ODM 


Primary subroutine which is called by the 
user's main program and supervises control 
of the iteration process. 

Initialize the default values of. control 
parameters and internal variables residing 
in the labeled COMMON blocks. 

Direction finding for unconstrained mini- 
mization process by means of modified 
Newton's method. If finite difference 
gradient calculation, is asked for, gradient 
of objective and constraint functions are 
computed in this subroutine. Also auto- 
mated search for an appropriate transition 
parameter is carried out whenever necessary. 
Called by ODM (one-dimensional minimization) 
and evaluates the composite function value 
for a given design. Quadratic extended 
penalty function scheme is administered by 
this subroutine. 

One-dimensional minimization is carried 
out along the direction provided by DIRCTN 
by means of the golden section algorithm. 
Linear constraints are treated separately 
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PRINTD 


RFACTR 


SAD007 


SAD008 


from nonlinear constraints, since the 
distances to linear constraint boundaries 
are readily calculated. 

This routine prints debugging information 
on a line printer. For users who wish to 
use NEWSUMT as a black box optimizer, this 
routine may be replaced by a dummary sub- 
routine which does nothing. 

Called only once at the beginning to compute 

the initial penalty multiplier r^. 

A part of the linear equation solver used 

in DIRCTN. A positive definite symmetric 

matrix stored in a vector form (upper tri- 

T 

angular only) is decomposed into a LDL 
form. If a matrix is almost singular or 
non-positive definite, an error flag is 
turned on, and the direction finding pro- 
cess is switched to the direction of steep- 
est descent. 

A part of the linear equation solver used 
in DIRCTN. Upon successful decomposition 
of the coefficient matrix of SADJ307 , this 
routine is called for back and forward 
substitutions , 
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V. 


PRACTICAL CONSIDERATIONS 


5 . 1 Formulation 

The standard form given in Eq. (1) is sufficient theore- 
tically, but for numerical stability and efficiency, it is 
always beneficial to normalize or to scale design variables and 
constraints. The optimization process tends to be stable when 
all Max(3g /3X.) have similar orders of magnitude and all 

- j q 3 

g (X) , q=l , 2 , . . . NTCE , also have similar orders of magnitude, 

4 

The user should be aware that ill conditioning of the Hessian 
matrix of <j> will result in extremely poor performance of this 
program. For example, consider a case where there is only one 
active constraint. Without contribution from the objective 
function, the rank of the matrix J of Eq. 9 is only 1. Haftka 
proposes to use diagonal perturbation (Ref. 2), but this also 
fails if the magnitudes of some diagonals are extremely small. 
Contributions from the objective function tend to alleviate this 
difficulty, provided they are significant, 

5 . 2 Choice of Control Parameters 

The most important parameter which the user should con- 
sider is the penalty multiplier decrease ratio (RACUT) . If 
RACUT is very small, then the penalty multiplier RA is decreased 
very rapidly and experience shows that iterative search tends 
to be terminated prematurely as the design gets very close to 
a small number of constraints. If RACUT is large, then conver- 
gence will be slow, but the optimal design will usually be of 
better quality. For difficult problems, RACUT can be as large 
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as 0 . 5~0 . 6 . For easy problems, RACUT may be taken as 0.1-0. 5. 
MAXRSF and RAMIN should be determined based on the assigned 
value for RACUT. 

The maximum allowable number of one-dimensional minimiza- 
tion per unconstrained minimization should be specified inde- 
pendent of the number of design variables. For most practical 
problems, the modified Newton's method is very effective and 
convergence will usually be achieved in three to five one- 
dimensional minimizations. The maximum number of golden section 
iterations (MAXGSN) is difficult to estimate, but the default 
value (20) has been used extensively with satisfactory perfor- 
mance. 

5 . 3 Variable Transfer Strategy 

The user may be uncomfortable to see a long list of argu- 
ments in SUBROUTINE NEWSUM. For ordinary small applications 
this looks awkward but when the problem size is large and 
optimum usage of main memory is required, this transfer strat- 
egy becomes beneficial since it allows dynamic array allocation. 
There are a number of cases where the analysis programs pro- 
vided by the user also require that certain arrays should be 
allocated dynamically. The arrays RAN and IAN are available 
for this purpose. Usually, all the arrays are packed in RAN 
and IAN sequentially and addressing pointers should be stored 
as a part of IAN. 
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VI. EXAMPLES 


Two simple examples are included to illustrate the scheme 
for writing computer programs to interface with NEWSUMT. They 
may also serve as test cases when NEWSUMT is initially imple- 
mented on a particular computer system. Both examples are self 
contained and need no input data, therefore an executable module 
created by linking each example program with the NEWSUMT pro- 
gram should readily be processed. Lineprinter outputs are 
abbreviated because of large volume of data from each iteration, 
but they should be sufficient to check the results which the 
user obtains. Even with single precision floating point vari- 
ables, on a 32 bits/word machine, all results should agree to 
at least three significant digits with the results given herein. 
In the printed output, TOTAL FUNCTION stands for the composite 
function <f>(X, r ) and OBJECTIVE FUNCTION refers to F (X) . 

IT 

(1) Minimize F(X) = lOX.^ + X 2 

Subject to g^X) = 2X 1 - X 2 - 1 2 0 

g 2 (X) = X 1 - 2X 2 + 1 A 0 
g 3 (X) = X^ 2 + 2X X + 2X 2 - 1 > 0 

This example is solved by a program given on p. 28 • The 
name of the subroutine for analysis is EXAMP 1 which must be 
declared in the main program as an external parameter. All 
gradient information is expressed analytically and computed in 
the analysis program as represented by vectors, DOBJ, DDOBJ 
and DG. 

An initial design for this particular example is chosen 
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as (2.0, 1.0), which satisfies all constraints and gives the 
objective function equal to 21.0. The final result given on 
page 20 indicates that an optimal design (0.5515, 0.1006) is 
obtained with corresponding objective function equal to 5.5917 
One can easily observe that critical constraints for this opti- 
mal design are g^ (X) and g^(X), since their values are 
0.3745x10 ^ and 0.4768x10 ^ while g 2 (X) is 1.351. 

This example is helpful in illustrating the iteration 
characteristics of the NEWSUMT optimizer. Figure 4 shows the 
design space and iteration paths using five different initial 
designs (one feasible, one critical, three infeasible). A 
conceptually attractive feature of the SUMT algorithm, namely 
its tendency to "funnel down the middle of the feasible space," 
is graphically illustrated by this example. It is this tendency 
of the NEWSUMT program, to generate a sequence of steadily 
improving noncritical designs, that make it especially well 
suited to engineering design applications. 

(2) Minimize F (X) = (X 1 2 -5X 1 ) + (X 2 2 -5X 2 ) + ( 2X 3 2 -21X 3 ) 

+(X 4 2 +7X 4 )+50 

g-L (x) = (-x 1 2 -x 1 ) + (-x 2 2 +x 2 ) + (-x 3 2 -x 3 ) 

+(-X 4 2 +X 4 )+8 - 0 
g 2 (^) = (-X 1 2 +X 1 ) + (-2X 2 2 ) + (-X 3 2 ) 

+(-2X 4 2 +X 4 )+10 - 0 
g 3 (X) = (-2X 1 2 -2X 1 ) + (-X 2 2 +X 2 ) + (-X 3 2 ) 

+X 4 +5“0 
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The optimal point of this well known Rosen-Suzuki problem 
is 

F (0. 0, 1.0, 2.0, -1.0) = 6.0 

The main program and the analysis program to solve this problem 
with NEWSUMT are given on pages 32 and 33, respectively. The 
analysis subroutine is called EXAMP 2 , which is declared as an 
external parameter in the main program. The subroutine EX AMP 2 
is written so that it can evaluate the analytic derivatives 
of objective and constraint functions. But in the particular 
run shown here, the capability to compute all gradients using 
finite difference is activated by specifying IDF = 5 in the 
main program, thus DOBJ, DDOBJ and DG are never computed in 
EXAMP2 . The user is encouraged to test the same program using 
analytic gradients, i.e. simply changing IDF to 0 in the main 
program. 

The initial design is selected as (1.0, 1.0, 1.0, 1.0) 
which satisfies all constraints and the corresponding objective 
function value is 31.0. The iteration process converges in 
10 stages, yielding an optimal design 

F ( 0 . 0020 , 1.000, 1.998, - 1.002) = 6.000427 

which agrees satisfactorily with the theoretical optimal design. 
Constraints g 1 (X) and g 2 (X) are active since they are respectively 
0.1545x10 ^ and 0.1106x10 ^ while g 2 (^) is 0,9983. 
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MAIN PROGRAM FOR PROBLEM EXAMPLE 1 


DECLARE THE NAME OF ANALYSIS PROGRAM AS 
EXTERNAL EXAMP1 


AN EXTERNAL 


ENTRY 


DIMENSION XO (2 ) , X(2) 

* DOBJC2) , DDOBJC3) 

* BU(2) , BL( 2 ) 

* G2(3) , G3C3) 

* FDCV(2) , RAN(l) 

COMMON/CONTRL/C , EPSGSN, 

* RA , RACUT , RAMIN , 

* IFD , JPRINT , JSIGNG, 

* MFLAG , NDV , NTCE 


, SNC2) , G( 3 ) 

» S( 2 ) , DH( 2 ) 

, ISIDEC 2 ) , GB( 3) 
, ILIN(3) , 

, IAN ( 1 ) 

EPSODM, EPSRSF, 
STEPMX , 

LOBJ , MAXGSN, 


, DG(6 ) , 

, IIK(2) , 
, G1 ( 3 ) 


GO , P 
MAXODM, MAXRSF 


100 


N1 = 2 
N2 = 3 
N3 = 3 
N4 = 6 
NDV = 2 
NTCE = 3 
LOBJ = 1 
DO 100 I = 
BL(I) = 0.0 
ILIN(I) = 1 
ISIDECI) = 
ILIJK3) = 0 
NRANDM = 1 
NIANDM = 1 
MFLAG = 0 
JPRINT = 1 



Objective, function it> LinaaA. 

Only lotoeA bound* o£ 0.0 a/i& impoizd 
to Xj and'xp. 

9 7 a nd g ^ aAe. tinoxvi, g ^ a/iz nonltmaJi. 



c 


c 


SPECIFY THE INITIAL DESIGN VARIABLES 
X0(1) = 2.0 
X 0 ( 2 ) = 1.0 


* 

* 

* 

* 

* 


CALL NEMSUMC EXAMP1 , 


BL 

t 

BU 

DDOBJ 

, DG 

DH 

, DOBJ 

f 


FDCV 

9 

FMIN , 

G 

, GB 

G1 

, G2 

9 G3 


OBJ 

9 

OBJMIN 9 

S 

, SN 

X 

, XO 

* 


IIK 

9 

ILIN , 

ISIDE 

, N1 

N2 

, N3 

, N4 


RAN 

9 

NRANDM, 

IAN 

, NIANDM : 

) 





STOP 

END 
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ANALYSIS PROGRAM FOR EXAMPLE PROBLEM 

SUBROUTINE EXAMPKINFO , X 

* GB , DG , N1 

* RAN , NRANDM, IAN 

DIMENSION X(N1 ) , DOBJ(Nl) , 

* DG(N4) , GB(N2) 


1 


OBJ , DOBJ , DDOBJ , 

N2 , N3 , N't 

NIANDM ) 

DD0BJ(N3 ) , G(N2) 

RAN (NRANDM ) , IAN ( NIANDM ) 


G 


t 


A = X(l) 
B = X(2) 


GO TO (100, 200, 300, 400, 500), INFO 

EVALUATE OBJECTIVE FUNCTION 
100 OBJ = 10.0*A + B 
RETURN 

EVALUATE CONSTRAINT FUCTIONS 
200 G ( 1 ) = 2.0*A - B - 1.0 
G ( 2 ) = A - 2.0*B + 1.0 
G( 3 ) = -A*A + 2 . 0*(A+B) - 1.0 
RETURN 

EVALUATE THE FIRST AND THE SECOND ORDER DERIVATIVES OF 

THE OBJECTIVE FUNCTION 

300 DOBJ ( 1 ) = 10.0 • ■ 

DOBJ ( 2 ) =1.0 
DDOBJ ( 1 ) = 0.0 
DDOBJ ( 2 ) = 0.0 
DDOBJ ( 3 ) = 0.0 
RETURN 

EVALUATE GRADIENT OF NONLINEAR CONSTRAINTS 
400 DG ( 3 ) = -2.0*A + 2.0 
DG ( 6 ) = 2.0 
RETURN 

EVALUATE GRADIENT OF LINEAR CONSTRAINTS 
500 DG ( 1 ) = 2.0 
DG( 2 ) = 1.0 
DG ( 4 ) = -1.0 
DG( 5 ) = -2.0 
RETURN 
END 


Note.: IF 0=0, t.e. analytical deAtvativeb fioa all fiunctionA 

one .computed . in the uaca'a analyidt paogam, ti 
the default option. 
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Example. PAoblem 1 - Abridged output [JPRINT=1] 


xxxxxxxxxxxxxxx******* NEWSUMT OPTIMIZER ****xx*xxx*xxxxx*xx*xx 


CONTROL PARAMETERS 

INITIAL TRANSITION POINT GO = 

TRANSITION POINT EXPONENT P = 

INITIAL TRANSITION POINT COEFFICIENT ... C = 

GOLDEN SECTION CONVERGENCE EPSGSN = 

UNCONSTRAINED MINIMIZATION CONVERGENCE EPSODM = 

CONVERGENCE AMONG RESPONSE SURFACES EPSRSF = 

RESPONSE FACTOR REDUCTION RATIO RACUT = 

MINIMUM ALLOWABLE RESPONSE FACTOR IRAMIN= 

MAXIMUM ALLOWABLESTEP SIZE STEPMX = 

MAXIMUM ALLOWABLE GOLDEN SECTIONS MAXGSN = 

MAXIMUM NUMBER OD O.D.M. PER SURFACE MAXODM = 

MAXIMUM ALLOWABLE RESPONSE SURFACES MAXRSF = 

PRINTOUT CONTROL JPRINT = 

FINITE DIFFERENCE GRADIENT CONTROL IFD = 


0.100 0E+ 0 0 
0 .5000E+00 
0.2000E+00 
0 . 1000E-02 
0. 1000E-02 
0 . 5000E-03 
0.1000E+00 
0 . 1000E-12 
0 . 1000E+11 
20 
6 


o Analytical. deAivativet 


SYSTEM PARAMETERS 


NUMBER OF DESIGN VARIABLES NDV = 2 

NUMBER OF EFFECTIVE CONSTRAINTS NTCE = 3 


INITIAL DESIGN ANALYSIS SUMMARY 

Xj and Xp 

LoiveA bound Aide conAtAaint £oa 
X p iA 1.000. 

0 . 2000 E + 01 0 . 1000E+01 O.lOOOE+Ol-i ^ ^ nc ^ ex ^ pOAitcoe, it 

objective function = o. 2100000 e+ 02 / tndtcatet uppeA bound Aide 

V / conAtAaint . ) 

"3j, 9 2 and 9 3 

- 

+JL IIAAJUK/L F(x”) 

) than itA value at the beginning 

ol tlvLA ACOAck. 

A 

P 

OPTIMIZATION OF RESPONSE SURFACE NO. 1 PENALTY MULTIPLIER = 0.525000E+01 



INITIAL DESIGN VARIABLE VECTOR 
0.2000E+01 0 . 1000E+01 

SIDE CONSTRAINTS 

-1 r =-2 

1 I 0 . 1 ( 


0 .2000E+0 
CONSTRAINTS 


lOOOE+Ol 




ONE DIMENSIONAL SEARCH 1 TOTAL FUNCTION= 0.399643E+02 

ONE DIMENSIONAL SEARCH 2 TOTAL FUNCTION= 0.398529E+02 

REPEAT ODM BY COMPLETE ANALYSIS 

ONE DIMENSIONAL SEARCH 3 TOTAL FUNCTION= 0.398509E+02 

ONE DIMENSIONAL SEARCH A TOTAL FUNCTION= 0.398347E+02 

RESULTS AT THE END OF THIS UNCONSTRaVeD MINIMIZATION 
INITIAL DESIGN VARIABLE VECTOR \ 

0 . 1557E+01 0 . 7895E+00 \ 

SIDE CONSTRAINTS \ , 

-1 -2 A ) 

0 . 1557E+0 1 0 . 7895E+00 ' 

CONSTRAINTS 

0 . 1322E+01 0 . 9786E+00 0.1269E+01 

TOTAL FUNCTION = 0 . 3983470E+02 

OBJECTIVE FUNCTION = 0 . 1633887E+02 


OBJECTIVE FUNCTION= 0.165580E+02 
OBJECTIVE FUNCTION= 0.162391E+02 


OBJECTIVE FUNCTION= 0.163712E+02 
OBJECTIVE FUNCTION= 0.163389E+02 
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OPTIMIZATION OF RESPONSE SURFACE NO. 2 


PENALTY MULTIPLIER = 0.525000E+00 


ONE DIMENSIONAL SEARCH 1 
ONE DIMENSIONAL SEARCH 2 


TOTAL FUNCTION 3 0.139056E+02 
TOTAL FUNCTION 3 0.136870E+02 


REPEAT ODM BY COMPLETE ANALYSIS 

ONE DIMENSIONAL SEARCH 3 TOTAL FUNCTION 3 0.136863E+02 
ONE DIMENSIONAL SEARCH 4 TOTAL FUNCTION 3 0:i36863E+02 

RESULTS AT THE END OF THIS UNCONSTRAINED MINIMIZATION 
INITIAL DESIGN VARIABLE VECTOR 
0 . 8574E+0 0 0 . 3626E + 00 

SIDE CONSTRAINTS 

-1 -2 
0.8574E+00 0 . 3626E + 00 

CONSTRAINTS 

0 . 3499E+0 0 0 . 1133E+01 0.7049E+00 

TOTAL FUNCTION = 0 . 1368631E+02 

OBJECTIVE FUNCTION 3 0 .8917752E+01 


OBJECTIVE FUNCTION 3 0.945919E+01 
OBJECTIVE FUNCTION 3 0.887688E+01 


OBJECTIVE FUNCTION 3 0.891434E+01 
OBJECTIVE FUNCTION 3 0.891775E+01 


OPTIMIZATION OF RESPONSE SURFACE NO. 10 PENALTY MULTIPLIER = 0.525000E-08 


ONE DIMENSIONAL SEARCH 1 
ONE DIMENSIONAL SEARCH 2 
ONE DIMENSIONAL SEARCH 3 


TOTAL FUNCTION 3 0.559644E+01 
TOTAL FUNCTION 3 0.559643E+01 
TOTAL FUNCTION 3 0.559642E+01 


OBJECTIVE FUNCTION 3 
OBJECTIVE FUNCTION 3 
OBJECTIVE FUNCTION 3 


RESULTS AT THE END OF THIS UNCONSTRAINED MINIMIZATION 
INITIAL DESIGN VARIABLE VECTOR J.UN 

0.5515E+00 0.1006E+00 

SIDE CONSTRAINTS 

-1 -2 
0.5515E+00 0 . 1006E+00 

CONSTRAINTS 

0 . 37 45E-04 0.1351E+01 0.4768E-04 

TOTAL FUNCTION = 0 . 5596422E+01 

OBJECTIVE FUNCTION 3 0 . 5596172E+01 


FINAL RESULTS OF OPTIMIZATION 


Optimal dteiqn 
= 0.5515 


1 

Fix 


3 0. 1006 
5.597672 


0 .559623E+01 
0 . 559616E+01 
0 .559617E+01 


CURRENT DESIGN VARIABLE VECTOR 
0 . 5515E+00 0 . 1006E + 00 

SIDE CONSTRAINTS 

-1 -2 
0.551 5E+ 0 0 0 . 1006E+00 

CONSTRAINTS 

0.3745E-04 0.1351E+01 0.4768E-04 

TOTAL FUNCTION 3 0 . 5596422E+ 0 1 

OBJECTIVE FUNCTION 3 0 . 5596 172E+0 1 


No 

3l 


6-ide. con^taalnti ajie cuttiaal. 
and g an. 2 . csiltlcal coiutaaliitb 
while, g ^ cum not. 


FINAL STATISTICS 


NUMBER OF RESPONSE SURFACE . . 
NUMBER OF ONE DIMENSIONAL SEARC 
NUMBER OF ANALYSES 


10 

34 


OBJECTIVE FUNCTION .... 

GRADIENT OF OBJECTIVE FUNCTION 

CONSTRAINT FUNCTIONS * | 

GRADIENT OF LINEAR CONSTRAINT FUNCTIONS* .* 
GRADIENT OF NONLINEAR CONSTRAINT FUNCTIONS 
APPROXIMATE CONSTRAINT FUNCTIONS 


309 

1 

110 

1 

34 

199 


1 1 

110+199 = 309 
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MAIM PROGRAM FOR PROBLEM EXAMPLE 2 


DECLARE THE NAME OF ANALYSIS PROGRAM AS EXTERNAL ENTRY 
EXTERNAL EXAMP2 


DIMENSION X0C4) 

* DGC12) , 

* DHU) 

* ISIDEU) , 

* G1C3) 

* RAN(l) 


X(4) 

DOBJU) 

IIKU) 

GB(3) 

G2(3) 

IAN(l) 


, SNC4) , G( 3 ) 

, DDOBJ(IO), S(4) 

, BUU) , BL(4) 

> ILINC3) , FDCVC4) 

, G3 ( 3 ) , XBC4) 


COMMON/CONTRL/ C , EPSGSN, EPSODM, EPSRSF, GO , P 

* * RA , RACUT , RAMIN , STEPMX 

* ' IFD » JPRINT , JSIGNG, LOBJ , MAXGSN, MAXODM, MAXRSF 

* * MFLAG , NDV , NTCE 


C 


C 


INITIALIZE NON-DEFAULT CONTROL PARAMETERS 
NRANDM = 1 
NIANDM = 1 
N1 = 

N2 = 3 

N3 = ( (Nl+1 ) * Nl)/2 
N't = N1 * N2 
NDV = 4 
NTCE = 3 
LOBJ = 0 
MFLAG = 0 
DO 100 I = 1, NDV 

100 ISIDE(I) = 0 

DO 110 I = 1, NTCE 

110 ILIN(I) =0 

IFD = 5 


JPRINT = 2 

INITIALIZE THE STARTING DESIGN 
DO 120 I = 1, NDV 
120 X0(I) = 1.0 


A tt vcuiiableA cme unbounded. 
(Wo t>tde con6tAcunt6) 


AtC comt/iainti and objective, 
function* cum nonltnean.. 

F twite dt^eAence scheme -Li, 
activated to compute. deAtvattve& 
o & alt nontineem function* . 


CALL NEWSUMC EXAMP2 


* 

, BL 

> 

BU 

DDOBJ 

9 

DG 

9 

DH 

, DOBJ 


* 

, FDCV 

9 

FMIN , 

G 

9 

GB 

9 

G1 

, G2 

, G3 

X 

, OBJ 

9 

OBJNIH, 

S 

9 

SN 

9 

X 

, XO 


* 

, UK 

9 

ILIM , 

ISIDE 

9 

HI 

9 

N2 

, N3 

, N4 

* 

, BAH 

9 

NRANDM, 

IAN 

9 

NIANDM) 





STOP 

END 
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ANALYSIS PROGRAM FOR EXAMPLE PROBLEM 2 


C 

c 


SUBROUTINE EXAMP2CINFO , X 
* g.„ ’ GB , DG 

RAN , NRANDM, IAN 


OBJ , 
N1 

NIANDM) 


DOBJ 

N2 


DDOBJ , 
N3 


N4 


DIMENSION X(4) 

* DG( 12) 

T = X(l) 

U = X(2) 

V = X ( 3 ) 

W = X(4) 


DOBJ (4 ) 
GB( 3 ) 


DDOBJ ( 10), G(3) , 

RAN(NRANDM) , IAN (NIANDM) 


C 

C 


C 

C 


C 

C 


C 

c 


GO TO (100, 200, 300, 400, 500), INFO 

2.0*V*V -21.0*V + W*W 


EVALUATE OBJECTIVE FUNCTION 
100 OBJ = T*T - 5 . 0*T + U*U - 5 0*U 
* + 7.0*W + 50.0 

RETURN 


C 

c 

c 


EVALUATE CONSTRAINT FUNCTIONS 
200 CONTINUE 

G ( 1 ) = -T*T - T - U*U + U - V*V 
G( 2 ) = -T*T + T - 2.0*U*U - V*V 
G ( 3 ) = -2.0*T*T - 2 . 0*T - U*U + 
RETURN 


- V - W*W 

- 2.0*W*W 
U - V*V + 


+ W + 8.0 
+ W + 10.0 
W + 5.0 


FUNCTION ™ E FIRST AND THE SEC0ND 
300 CONTINUE 

DOBJ ( 1 ) = 2.0*T 
DOBJ ( 2 ) = 2.0*U 
DOBJ ( 3 ) = 4.0*V 
DOBJ ( 4 ) = 2.0*W 


5.0 

5.0 

21.0 

7.0 


310 


DO 310 I 
DDOBJ (I ) 
DDOBJ ( 1 ) 
DDOBJ ( 3 ) 
DDOBJ ( 6 ) 
DDOBJ (10) = 
RETURN 


1, 10 
0 . 0 
2.0 
2.0 
4 . 0 
2.0 


ORDER DERIVATIVES OF OBJECTIVE 


Note.: TlvLi, pant ol the. pnognam 

li> not uted, &lnee IFP=5 
-ci -ipectfited tn the MAIW 

pnognam. 


EVALUATE GRADIENT OF NONLINEAR Coh 
400 CONTINUE 


DG ( 1 ) 
DG ( 2 ) 
DG( 3 ) 
DG ( 4 ) 
DG( 5 ) 
DG ( 6 ) 
DG ( 7 ) 
DG ( 8 ) 
DG( 9 ) 


-2 . 0*T 
-2. 0*T 
-4. 0*T 
-2. 0*U 
-4 . 0*U 
_ 2 . 0*U 
-2. 0*V 
-2. 0*V 
-2 . 0*V 


1.0 

1.0 

2.0 

1.0 

1.0 

1.0 


DG(10) = -2.0*M 
DG(ll) = -4 . 0*W 
DG( 12 ) = 1.0 
RETURN 


1.0 

1.0 


ISTRAINTS 


5oo A continue ADIENT ° F LINEAR constraints 

RETURN 
END 
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Example. Problem 2 - Abridged Output ( J PRINT =2 ) 


NEWSUMT OPTIMIZER 


******XXX*XXXXXXXXXXXX 


CONTROL PARAMETERS 

INITIAL TRANSITION POINT GO = 

TRANSITION POINT EXPONENT P = 

INITIAL TRANSITION POINT COEFFICIENT ... C = 

GOLDEN SECTION CONVERGENCE EPSGSN = 

UNCONSTRAINED MINIMIZATION CONVERGENCE EPSODM = 

CONVERGENCE AMONG RESPONSE SURFACES EPSRSF = 

RESPONSE FACTOR REDUCTION RATIO RACUT 3 

MINIMUM ALLOWABLE RESPONSE FACTOR IRAMIN= 

MAXIMUM ALLOWABLESTEP SIZE STEPMX = 

MAXIMUM ALLOWABLE GOLDEN SECTIONS MAXGSN = 

MAXIMUM NUMBER OD O.D.M. PER SURFACE MAXODM = 

MAXIMUM ALLOWABLE RESPONSE SURFACES MAXRSF = 

PRINTOUT CONTROL JPRINT = 

FINITE DIFFERENCE GRADIENT CONTROL IFD = 


0.1000E+00 
0 . 5000E+00 
0 . 2000E+00 
0 . 1000E-02 
0.1000E-02 
0.5000E-03 
0 . 1000E+00 
0 . 1000E-12 
0.100 0E+1 1 
20 
6 
30 
2 
5 


SYSTEM PARAMETERS 


NUMBER OF DESIGN VARIABLES NDV = 4 

NUMBER OF EFFECTIVE CONSTRAINTS NTCE = 3 


INITIAL DESIGN ANALYSIS SUMMARY 

INITIAL DESIGN VARIABLE VECTOR 

0.1000E+01 0 . 1000E+01 0.1000E+01 0.1000E+01 

CONSTRAINTS 

0.4000E+01 0.6000E+01 0.1000E+01 

OBJECTIVE FUNCTION = 0 . 31 0 0 0 0 0E+02 


OPTIMIZATION OF RESPONSE SURFACE NO. 1 PENALTY MULTIPLIER = 0.218824E+02 


DIRECTION FINDING 

TRANSITION POINT = 0.100000E+00 


^ECTION COMPUTED BY MODIFIED NEWTON'S METHOD 

oLU t h - - 0 . 226 4E+ 0 2 

NORMALIZED DIRECTION VECTOR 
-0.2707E+00 -0.2924E+00 0.3743E+00 -0.8374E+00 


-- ONE DIMENSIONAL MINIMIZATION RUN NO 1 
-- END OF O.D.M. DISTANCE FOR MIN. PT. = ( 
ONE DIMENSIONAL SEARCH 1 TOTAL FUNCTION 3 


^ 

9407793E+00 
. 513895E+02 



OBJECTIVE FUNCTION= 


0.205227E+02 


NOT CONVERGED - CHECK1 3 -0.1000E+01 CHECK2 3 0.1946E+19 EPSODM= 0.1000E-02 

DIRECTION FINDING 

TRANSITION POINT = 0 . 425813E+00 


SECTION COMPUTED BY MODIFIED NEWTON'S METHOD 
SLOPE = -0.1841E+02 
NORMALIZED DIRECTION VECTOR 
-0.3202E+00 0.9768E-01 0.3788E-01 -0.9416E+00 


-- ONE DIMENSIONAL MINIMIZATION RUN NO. 2 
— END OF O.D.M. DISTANCE FOR MIN. PT. = 0 
ONE DIMENSIONAL SEARCH 2 TOTAL FUNCTION 3 


.8987486E+00 

0.434400E+02 


OBJECTIVE FUNCTION 3 0.152161E+02 
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OPTIMIZATION OF RESPONSE SURFACE NO. 10 


PENALTY MULTIPLIER = 0.218824E-07 


- DIRECTION FINDING 

TRANSITION POINT = 0.912702E-05 

DIRECTION COMPUTED BY MODIFIED NEWTON'S METHOD 
SLOPE = -0 . 9156E+01 
NORMALIZED DIRECTION VECTOR 

0 . 4763E+00 0 . 7524E+00 0.3046E+00 -0.3379E+00 

- ONE DIMENSIONAL MINIMIZATION RUN NO. 1 

- END OF O.D.M. DISTANCE FOR MIN. PT. = 0 . 6073422E-04 

ONE DIMENSIONAL SEARCH 1 TOTAL FUNCTION= 0. 60007 9E+01 OBJECTIVE FUNCTION= 0.600052E+01 


ONE DIMENSIONAL MINIMIZATION RUN NO. 3 


REPEAT ODM BY COMPLETE ANALYSIS 

- END OF O.D.M. DISTANCE FOR MIN. PT = 0 1 7 r 

ONE DIMENSIONAL SEARCH 3 TOTAL FUNCTION= 6 . 600077E+01 


OBJECTIVE FUNCTION= 0 . 6 0 0 043E+0 1 


CONVERGED CHECI<1= 0.2861E-05 CHECK2= 0.9537E-06 EPSODM= 


RESULTS AT THE END OF THIS UNCONSTRAINED MINIMIZATION 
INITIAL DESIGN VARIABLE VECTOR 

CONSTRAINTS 2 °- 1000E+01 0-1998E+01 -0.1002E+01 

0-1545E-03 0 . 9983E+00 0.1106E-03 

FUMCTI0N = 0.6000766E+01 

OBJECTIVE FUNCTION = 0 . 6000427E+01 

(F(X ' 


0.1000E-02 


10 


- F(X 9 )) / FIXg) 


CONVERGED - CHECK3= 0.4368E-03 CHECK4= 

DOUBLE CONVERGENCE CRITERIA IS SATISFIED 
FINAL RESULTS OF OPTIMIZATION 


CURRENT DESIGN VARIABLE VECTOR 

CONSTRAINTS 2 °- 1000E+01 °-1998E+01 -0.1002E+01 

0 . 1545E-0 3 0 . 9983E+00 0.1106E-03 

TOTAL FUNCTION = 0 . 60 00766E+0 1 

OBJECTIVE FUNCTION = 0 . 6 00 0427E+0 1 

FINAL STATISTICS 

NUMBER OF RESPONSE SURFACE .... 10 

NUMBER OF ONE DIMENSIONAL SEARCH . 35 

NUMBER OF ANALYSES 

OBJECTIVE FUNCTION .... 

GRADIENT OF OBJECTIVE FUNCTION 

CONSTRAINT FUNCTIONS . ' ' 

GRADIENT OF LINEAR CONSTRAINT FUNCTIONS' ! 
GRADIENT OF NONLINEAR CONSTRAINT FUNCTIONS 
APPROXIMATE CONSTRAINT FUNCTIONS .... 


0.1163E-03 EPSRSF= 0.5000E-03 

(F (X J0 ) - F(X 9 )J / F(X 9 ) 

Both CHECK3 and CHECK4 oj le irr.Mvi 
than EP5RSF. 


699 

0 

285 

0 

0 

274 


35 



REFERENCES 


1. Fiacco, A.V. and McCormick, G.P., "Nonlinear Programming: 
Sequential Unconstrained Minimization Techniques," 

John Wiley and Sons, Inc., 1968. 

2. Haftka, R. , "Automated Procedure for Design of Wing Struc- 
tures to Satisfy Strength and Flutter Requirements," 

NASA TN D-7264 , 1973. 

3. Schmit, L.A. and Miura, H. , "An Advanced Structural 
Analysis/Synthesis Capability - ACCESS 2," Int. J. Num. 
Meth. Engrg . , Vol. 12, No. 2, February 1978, pp. 353-377. 

4. Kavlie, D. and Moe, J., "Automated Design of Frame Struc- 
tures," J. of the Structural Div., ASCE, Vol. 97, No. ST1 , 
January 1971, pp. 33-62. 

5. Haftka, R. and Starnes, J.H., "Applications of a Quadratic 
Extended Interior Penalty Function for Structural Optimiza- 
tion," AIAA Journal, Vol. 14, No. 6, June 1976, pp. 718-724 

6. Cassis, J. and Schmit, L.A. , "On Implementation of the 
Extended Interior Penalty Function," Int. J. Num. Meth. 
Engrg., Vol. 10, No. 1, January 1976,- pp. 3-23. 


36 




* Only this subroutine is called by the user's 
main program to activate the NEWSUMT program. 


Fig. 1 Basic Program Organization 
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Fig. 2 Structure of the Analysis Program to Evaluate 
All Functions and Their Derivatives (if available) 
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Initial Design 

End of Unconstrained Minimization 




Optimal Design 


Fig. 4 Iteration Trajectories for Example 1. 
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Golden section algorithm is initiated using 
the function values evaluated at a 2 , a 4 , 

Fig. 5 One Dimensional Search Scheme 





APPENDIX A 


An Interpretation of Quadratic Extended Penalty Functions 

The NEWSUMT program is written based on the philosophy of 
interior penalty functions. However, implementation of the 
quadratic extended penalty function has made it possible to 
regard this program as an integration of interior and exterior 
penalty function philosophy. Let the composite function of an 
exterior penalty function be defined as 


where 


+ -v Q 

<f> Fy (X,R ) = F (X) + R V K (X) 

EX p p ^ q 

( 0 if g (x)> o 

v $) = ! 

( [g q (S>] 2 if V x)< o 


(Al) 


(A2) 


Now compare this with a simplified version of Eqs. (3) and (4) 


*IN ( . X ' r p> " 


F (X) + r 


Q 

l 

q=i 


H (X) 

q 


(A3) 



( A4 ) 


If a constraint is violated significantly, the dominant term in 
H (X) is clearly the first term 

H. 

H q (X)c[g g (X) ] 2 /e 3 (A5) 
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Recalling Eq. (7) for the relation e=Cv / r 
" P 


H (X) s [g (X)] 2 /C 3 r 3/2 (A6) 

Si p 

Therefore, if R p in Eq. (Al) is chosen as 

R p = l/(C 3 r p 1/2 ) (A7 ) 

then the relation given by Eq. (A4) may be interpreted so that 
it includes the basic ingredients of Eq. (A2) . This is illus- 
trated m Fig. 3. it is well recognized that both interior and 
exterior penalty functions exhibit poor numerical behaviors at 
or near the constraint boundaries where g q (X)=0. The quadratic 
extended penalty function interpolates the penalty function 
adequately in this critical region. 
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APPENDIX B 


CPU Timing Routine 

The function CTIME(I) is used to measure CPU time spent 
for various parts of data processing in the NEWSUMT program. 
Timing routines are usually installation dependent and CTIME(I) 
is available only on IBM 360/91 at UCLA. The function CTIME(I) 
gives the remaining CPU time, which is the difference between 
the estimated maximum CPU time specified in the JOB card and 
the CPU time already spent on the particular job. 

For most CDC computers, a function SECOND does a similar 
job, but it gives the CPU time expended by the particular job. 
Therefore, CDC users may add a simple function such as 

FUNCTION CTIME(I) 

CALL SECOND (T) 

CTIME=-T 

RETURN 

END. 
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